Behavior of the S Parameter in the Crossover Region Between Walking and 
QCD-Like Regimes of an SU(A^) Gauge Theory 
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We consider a vectorial, confining SU(A'^) gauge theory with a variable number, Nf, of massless 
fermions transforming according to the fundamental representation. Using the Schwinger-Dyson and 
Bethe-Salpeter equations, we calculate the S parameter in terms of the current-current correlation 
functions. We focus on values of Nf such that the theory is in the crossover region between the 
regimes of walking behavior and QCD-like (non-walking) behavior. Our calculations indicate that 
the contribution to S from a given fermion decreases as one moves from the QCD-like to the walking 
regimes. The implications of this result for technicolor theories are discussed. 

PACS numbers: PACS No. 



I. INTRODUCTION 

The properties of a vectorial gauge theory as a function 
of the fermion content are of fundamental importance. 
Here we consider such a theory (in (3-1-1) dimensions 
at zero temperature and chemical potential) with gauge 
group SU(7V) and Nf massless fermions transforming ac- 
cording to the fundamental representation of this group. 
For N — Nc — 3 and Nf = 2, this is an approximation to 
actual quantum chromodynamics (QCD) with just the u 
and d quarks, since their current-quark masses are much 
smaller than the scale Aqcd — 400 MeV. We restrict 
here to the range Nf < (11/2)A^ for which the theory 
is asymptotically free. An analysis using the two-loop 
beta function and Schwinger-Dyson equation (reviewed 
below) leads to the inference that for Nf in this range, 
the theory includes two phases: (i) for < Nf < Nf^cr a 
phase with confinement and spontaneous chiral symme- 
try breaking (SxSB), and (ii) for 7V/,„- < ^/ < (ll/2)iV 
a non-Abelian Coulomb phase with no confinement or 
spontaneous chiral symmetry breaking. We shall refer to 
Nf,cr, the critical value of Nf, as the boundary of the 
non-Abelian Coulomb (conformal) phase 1]. 

For Nf slightly less than Nf^cr, the theory exhibits an 
approximate infrared (IR) fixed point. Let the SU(A^) 
running gauge coupling be denoted as g{n), where ^ de- 
notes the energy or momentum scale, and let a(/i) = 
^(/i)^/ (47r). As /X decreases from large values, a{p) grows 
to be 0(1) at a scale A, but increases only rather slowly as 
fi decreases below A, so that there is an extended interval 
in energy below A where a is large, but slowly running 
("walking"). In addition to its intrinsic field-theoretic 
interest, this walking behavior is an essential ingredient 
of modern technicolor models of dynamical electroweak 
symmetry breaking 0|, providing the requisite enhance- 
ment of Standard-Model (SM) fermion masses As 
Nf approaches Nf^cr from below, quantities with dimen- 
sions of mass vanish continuously; i.e., the chiral phase 
transition separating phases (i) and (ii) is continuous. 

In this paper, we shall use solutions of the Schwinger- 



Dyson (SD) and Bethe-Salpeter (BS) equations to com- 
pute a derivative of the difference of the vector and axial- 
vector current-current correlation functions, nyy(O) — 
n^^(O). Up to a multiplicative factor, this is the coef- 
ficient Lio of one of the terms in the effective chiral La- 
grangian for the theory |l[Tol|. Moreover, in the context 
in which one considers the SU(A^) theory as a technicolor 
(TC) model, with N ~ Ntc, the above quantity is pro- 
portional to the correction to the Z propagator due to 
virtual electroweak-nonsinglet technicolor particles, of- 
ten denoted as 5 0, E II3- We focus on the 
crossover region between the walking regime that occurs 
for Nf ^ Nf^cr and the QCD-like (non-walking) regime 
that occurs for smaller Nf. 

There are several motivations for this work. The quan- 
tity S, or equivalently Lio, is an intrinsic property of the 
SU(A^) theory, and it is of interest to understand how 
this quantity depends on Nf. Further, our calculations 
have important implications for technicolor models of dy- 
namical electroweak symmetry breaking Q-Q- For this 
application, as noted, we identify the SU(A'^) group with 
the technicolor gauge group. Precision electroweak data 
[TgI \l% determine allowed regions of values of S and the 
other Z and W propagator corrections denoted T and 
U, and yield a stringent constraint on the contributions 
from new particles in technicolor models. In order to as- 
sess the viability of these models, it is necessary to have a 
reliable calculation of the contribution to S from techni- 
color particles. A perturbative calculation of S is not re- 
liable since the technifermions are strongly interacting on 
the relevant scale ~ mz. Although it is possible to carry 
out a nonperturbative calculation for technicolor theories 
that behave like scaled-up QCD (by using QCD data as 
input for the relevant spectral functions), such theories 
are excluded since they cannot produce sufficiently large 
standard-model fermion masses. It is more difficult to 
carry out a nonperturbative estimate of S for technicolor 
models that have walking behavior; such estimates have 
been presented in Refs. 19J- 23] . These suggest that the 
contribution to 5' per technifermion electroweak doublet 
is reduced in the walking region. Meson masses and fp 



(the generalization of /^) were calculated in this walking 
regime in Ref. (See also Ref. [22| for the analogous 

calculations in QCD.) A motivation for our calculations 
in this paper is to gain further insight into the behav- 
ior of S by studying its behavior in the crossover region 
between the walking and QCD-like regimes. One of the 
reasons for concentrating on this crossover region is that 
the theory still has an approximate infrared fixed point 
here, and hence one can study the dependence of S on Nf 
without having to introduce a model-dependent cutoff on 
the growth of the SU(Af) gauge coupling in the infrared 
that was required in calculations of S via the SD and BS 
equations for small Nf values, such as Nf = 3 in QCD 
|2ll . As our previous calculations of meson masses 
and fp in this crossover region showed, [2^, although it 
is a restricted interval in Nf or equivalently, the value of 
the infrared fixed point, it is sufhcient large to observe a 
significant change between walking and QCD-like behav- 
ior. 

This paper is organized as follows. Section II is de- 
voted to a review some background material concern- 
ing the beta function, approximate infrared fixed point 
and walking behavior, and technicolor models. Section 
III contains definitions of the current-current correlation 
functions and related spectral functions together with the 
expression for S in terms of these correlation functions 
and an equivalent integral of spectral functions. In this 
section we also give the definition of Liq as a coefficient 
of a certain operator in the low-energy effective chiral 
Lagrangian for the theory. Section IV explains how the 
current correlators are obtained from Bethe-Salpeter am- 
plitudes while Section V discusses the method of solution 
of the Bethe-Salpeter equation. Our results are presented 
in Section VI and their implications for technicolor theo- 
ries in Section VII. An appendix discusses some analytic 
approximations. 



II. SOME PRELIMINARIES 

In this section we review some background for our cal- 
culations. For the theory under consideration, with an 
SU(iV) gauge group and Nf massless fermions in the 
fundamental representation, the renormalization group 
(RG) equation for the running coupling a(/i) is 



f3 = ti 



da{fi) oi{ixY 
dfj, 27r 



An 



(2.1) 



where /i is the momentum scale. The two terms listed are 
scheme- independent. The next two higher-order terms 
have also been calculated but are scheme-dependent; 
their inclusion does not significantly affect our results. 
For the relevant case of an asymptotically free theory, 
&o > so that an infrared fixed point exists if and only if 
bi < 0. This coefhcient bi is positive for < Nf < Nfjn, 
where Nfjp = (34iV^)/(13iV^ - 3), and negative for 
larger Nf. For TV = 3, Nfjn ~ 8.1 H^. The value 
of a at this IR fixed point, denoted a*, is given by 



a* = — 47rbn/&i ■ Substituting the known values of these 
terms ISO, Isil , one has 



-47r(llA^- 2A^/) 



34iV2 - ISNNf + 3N~^Nf ' 
Solving Eq. H2.2|l for Nf in terms of a* yields 
2N^[l7N{a^/TT) + 22] 



N 



f 



{13N^-3)(a^/TT) + 8N 



(2.2) 



(2.3) 



As is evident from Eqs. H2.2|l and (|2.3|l . a* is a mono- 
tonically decreasing function of Nf and equivalently Nf 
is a monotonically decreasing function of a*, for Nfjp, < 
Nf < (ll/2)iV. 

To study the dependence of S on Nf, what we actually 
vary is the value of the approximate IR fixed point a*, 
which depends parametrically on Nf. For definiteness, 
we shall take iV = 3; however, as will be seen, N only 
enters indirectly, via the dependence of the value of the 
infrared fixed point a* (Eq. (|2.6I) below) on N^.. Hence, 
our findings may also be applied in a straightforward way, 
with appropriate changes in the value of a*, to an SU(iV) 
gauge theory with a different value of N. 

In the one-gluon exchange approximation, the 
Schwinger-Dyson gap equation for the inverse propagator 
of a fermion transforming according to the representation 
R of SU(iV) has a nonzero solution for the dynamically 
generated fermion mass, which is an order parameter for 
spontaneous chiral symmetry breaking, if a > acr, where 
acr is given by 



1, 



(2.4) 



and C2(i?) denotes the quadratic Casimir invariant for 
the representation R 1^3 ■ Using 



C2{fund.) = C2f 



N^ 



1 



2N 



for the fundamental representation yields 

27rA^ 



3(7V2 _ 1) ■ 



(2.5) 



(2.6) 



For the case = 3 that we use for definiteness here, Eq. 
(|2.t)|) gives acr = t^I'^ — 0.79. To estimate A^/^.r, one 
solves the equation a* = acr, yielding the result 



A^. 



f-cr 



2N{50N^ - 33) 
5(5A^2 _ 3) 



(2.7) 



For the values A^ = 3 and N ^ 2 this gives A^/,cr — 
11.9 and Nf^cr — 7.9, respectively. These estimates are 
only rough, in view of the strongly coupled nature of the 
physics. Effects of higher-order gluon exchanges have 
been studied in Ref. [33 ■ In principle, lattice gauge 
simulations should provide a way to determine A^/,crj but 
the groups that have studied this have not reached a 
consensus |3J|. 



2 



If a* is less than the critical value, ttcr, for a bilinear 
fermion condensate to form, the above IR fixed point is 
exact, with the coupling a approaching a* from below as 
the energy scale E decreases from large values to zero. 
Let us denote the fermions as /f with a = 1, iV and i — 
1, ...,Nf. If > acr, as E descends from large values, 
the coupling a eventually exceeds the above critical value, 
the fermion condensates {J2a=i fa.ifi) (with no sum on 
i) form, and are equal for each flavor i — 1, Nf (with 
electroweak interactions negligibly small relative to the 
SU(iV) interaction). Accordingly, the global SU(A^/)i x 
SlJ{Nf)ji X U(l)y symmetry (where U(l)y is fermion 
number) is broken to generalized isospin times fermion 
number, SU(7V/)y x f7(l)y. 

Associated with this, the fermions pick up dynamical 
masses E and are integrated out as the energy scale de- 
creases below S. Hence, in this case, the IR fixed point 
is only approximate since in the effective field theory for 
energies below E, the form of the beta function is that 
for the pure gauge theory, Nf — 0. The case where a, 
is close to, and slightly larger than, acr, yields walk- 
ing behavior. In the strong walking regime, the dynam- 
ical fermion mass, and also hadron masses are exponen- 
tially smaller than the scale A at which the coupling first 
becomes 0(1) as the energy scale decreases from large 
values. Although our SD and BS equations are semi- 
perturbative, the analysis is self-consistent in the sense 
that our acr really is the value at which, in our approx- 
imation, one passes from the confinement phase to the 
non-Abelian Coulomb phase, and our values of a do span 
the interval over which there is a crossover from walking 
to QCD-like (i.e., non- walking) behavior. 

As is evident from the above results, decreasing Nf 
below Nf^cr has the effect of increasing a* and thus mov- 
ing the theory deeper in the phase with confinement and 
spontaneous chiral symmetry breaking, away from the 
boundary with the non-Abelian Coulomb phase. This is 
the key parametric dependence that we shall use for our 
study. In Refs. f24>.26j the range of a, used for the calcu- 
lation of meson masses was chosen to be 0.89 < a, < 1.0, 
an interval where there is pronounced walking behavior. 
For the case = 3 considered in Ref. 24] and here, 
given the above-mentioned value, acr — 7r/4, it follows 
that this lower hmit, a* = 0.89, is about 12 % greater 
than this critical coupling. The reason for this choice of 
lower limit on a, was that the calculation of S involves 
very strong cancellations as a* — acr 0+, rendering 
it progressively more and more difficult to obtain accu- 
rate numerical results in this extreme walking limit. For 
our present study of S we consider an interval extend- 
ing to larger couplings, from a* = 1.0 to a* = 1.8. Our 
upper limit is chosen in order for the ladder approxima- 
tion used in our solutions of the Schwinger-Dyson and 
Bethe-Salpeter equations to have reasonable reliability. 
From Eq. H2.3() it follows that a, — 0.89 corresponds to 
Nf = 11.65, about 2 % less than Nf^cr- For a couphng 
as large as a* = 1.8, the semi-perturbative methods used 
to derive eqs. (|2.2|l and (|2.3|l are subject to large cor- 



rections from higher-order perturbative, and from non- 
perturbative, contributions; recognizing this, the above 
upper limit of a* corresponds to Nf ~ 10.3, a roughly 13 
% reduction from Nf^cr = 11.9. Although this shift ap- 
pears to be by only a modest amount when expressed in 
terms of Nf^ in terms of a* it is a factor of two, and our 
calculations in Ref. 25] showed a dramatic change in the 
values of meson mass ratios and fp/A in this range, with 
these values changing from their walking limits toward 
QCD-like values. Hence we anticipate that this range 
can be sufficient to study the shift in the value of 5*, and 
our results confirm this. 

If Nf > Nf^cr, i-S., a* < Ucr so that this IR fixed point 
of the two- loop RG equation is exact, then, denoting 6 = 
6o/(27r), the solution to this equation can be explicitly 
written |35t . .36] in the entire energy region as 



a(At) 



w{e-\^,/kf 



1 



(2.8) 



where W{x) — F ^{x), with F{x) = xe^ , is the Lambert 
W function, and A is a RG-invariant scale defined by 



A 



/i exp 



1 



■In 



a* 



a{fi) 



a{fi) 



1 



a(/i) 



(2-9) 

Now since we are studying the confined phase with 
Nf < Nf^cn (q^* > cicr) with spontaneous chiral sym- 
metry breaking, a* is only an approximate, rather than 
exact, IR fixed point. Hence, the solution (|2.8(l is only 
applicable in an approximate manner to our case; for 
momenta much less than the dynamical fermion mass E, 
the fermions decouple, and in this very low-momentum 
region, with the fermions integrated out, the resultant a 
would increase above the value a* at the approximate 
IR fixed point. However, since E <C A in a walking or 
near-walking theory, it follows that this lowest range of 
momenta makes a small contribution to the relevant in- 
tegrals to be evaluated in our calculations. Hence, over 
most of the integration range for these integrals where 
the coupling a is large, it is approximately constant and 
equal to its fixed-point value, a* (see Fig. 2 of Ref. 26]). 
This means that one can use, as a reasonable approxima- 
tion, the expression 



a{^) — a* 9{A — ii), 



(2.10) 



where 9 is the step function. (This is the same approxi- 
mation used in Refs. H^H^.) Thus, in the walking re- 
gion and the adjacent crossover region, the calculations 
have the advantage that one can avoid having to intro- 
duce an artificial cutoff on the growth of a in the infrared, 
in contrast to the situation for smaller Nf, where the 
walking behavior disappears and this cutoff is necessary. 

Since an important application of our results is to tech- 
nicolor models, we briefly mention some relevant fea- 
tures of these models. As noted, the technicolor gauge 
theory has a gauge group SU(A^tc) and an asysmp- 
totically free coupling that gets large at the TeV scale 
13]. It contains a set of massless, vectorially coupled 
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technifermions. The left-handed components of these 
fermions transform as doublets under SU(2)i. The spon- 
taneous chiral symmetry breaking and formation of a bi- 
linear technifermion condensate breaks the electroweak 
symmetry from SU(2)l x U(l)y to U(l)em, produc- 
ing masses for the W and Z given to leading order by 
= m| cos^ 0w = (g'^ /A)Nd,tfF^(j, where Ftc is 
the technicolor analogue of /^r. In order to give masses 
to the Standard-Model fermions (which are technisin- 
glets), it is necessary to embed technicolor in a larger, 
extended technicolor (ETC) theory 113, 113, with inter- 
actions that transform technifermions to the Standard- 
Model fermions and vice versa. To satisfy constraints 
from flavor-changing neutral-current (FCNC) processes, 
the ETC vector bosons that mediate generation-changing 
transitions must have large masses, ranging from a few 
TeV to 10'^ TeV. For our present study, concerned with 
5, we concentrate on the technicolor theory at the scale 
of a few hundred GeV, with the ETC gauge bosons inte- 
grated out. 

We focus here on models in which the technifermions 
transform according to the fundamental representation of 
the S\]{Ntc) gauge group. Two simple examples are the 
so-called one-doublet and one-family technicolor models. 
A one-doublet TC model has Nf — 2 technifermions, 
denoted U and D, whose chiral components transform 
according to 



(iV,l,2). 



1,L 



Fl 
Ur 



: (iVTc,l,2)o,L, 



Dr : (iVTC,l,l)- 



l.R 



(2.11) 



where the numbers in parentheses refer to the dimensions 
of the representations of SU(A^tc) x SU(3)c x SU(2) l and 
the subscripts refer to the weak U(l)y hypercharges. The 
value Ntc = 2 has been preferred in recent TC/ETC 
model-building 0,E2 for several reasons, including the 
fact that it (i) minimizes technicolor contributions to the 
S parameter, (ii) can naturally produce a walking the- 
ory in a one- family model (see below), and (iii) makes 
possible a mechanism to explain light neutrino masses 
li^. The one-doublet TC model has one SU(2)l dou- 
blet of technifermions for each technicolor index, which 
we express as Nd^tf = 1, and hence a total number 
of SU(2)i doublets of technifermions equal to No^tot — 
Nd,tfNtc = Ntc- The TC sector with just these 
Nf = 2 technifermions would not exhibit walking behav- 
ior, but one can add SM-singlet technifermions to pro- 
duce a theory that does have such behavior jl^. In a 
one-family TC model the technifermions transform as 

Ql: (A^,3,2)i/3,i 
UR : {N,3,1)4/3^R 

dR-. (^,3,l)_2/3,i? 



Nr: (iV,l,l)o,fl 
Er: {N,1,1)-2J 



(2.12) 



Hence, this type of technicolor models contains Nf — 
2{Nc -I- 1) = 8 technifermions. As is evident from Eq. 

with N = Ntc = 2, the value iV/ = 8 is close to the 
value Nf^cr and hence, to within the accuracy of the two- 
loop beta function analysis, this technicolor model can 
naturally exhibit walking behavior. Reverting to general 
N = Ntc for our discussion, this one-family technicolor 
model thus has Nd,tf = (^c + 1) = 4 SU(2)l doublets 
for each technicolor index, and hence a total of Nrj^tot. = 
ANtc SU(2)l doublets of technifermions. 

Because of the spontaneous chiral symmetry breaking 
in the technicolor theory, the technifermions pick up dy- 

— 1/2 

namical masses Stc proportional to Ftc^tc ' where 
we have included the A'T^c'-dependent factor that would 
be present in the large- A'rc limit, since fp and Stc 
scale, respectively, like N^^ and Nj,q in this limit. For 
the one-doublet and one-family TC models, Ftc — 250 
and 125 GeV, respectively. In QCD, the constituent 
quark mass S ~ 3.5/^, and one expects a roughly simi- 
lar ratio in TC theories (see Fig. 3 in our previous work 
|25|). Since the SM gauge couphngs are small at the 
technicolor scale, different technifermions are expected to 
have roughly degenerate dynamical masses, and the con- 
tributions of the techniquark and technilepton doublets 
to one-loop corrections to the Z propagator are approxi- 
mately equal. 



III. EXPRESSION FOR S IN TERMS OF 
CURRENT-CURRENT CORRELATION 
FUNCTIONS 

As a measure of corrections to the Z propagator aris- 
ing from heavy particles in theories beyond the Standard 
Model, S was originally defined as [lj| 



4s 



dlifPi,-) 



i{mz) 



dq^ 



(3.1) 



s2 



1 



"W 



sin dvv, evaluated at mz and 
the superscript NP refers to the fact that the definition 
includes new physics beyond the Standard Model. In the 
case of technicolor, the technifermions are taken to have 
zero masses; because of the spontaneous chiral symme- 
try breaking in the TC theory, they pick up dynamical 
masses Stc of order the technicolor scale. More recent 
analyses of precision electroweak data define S slightly 
differently, replacing the derivative at = by a dis- 
crete difference (in the MS scheme) [ig 



S 



As 



PDG 



limz) 



n 



'zT(ml)- 



n 



(NP) 
ZZ 



(0) 



(3.2) 
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The difference between these definitions is small if the 
heavy physics scale Etc satisfies {2Titc l^zY ^ ^^ '"^ 
is the case in the TC models considered here. 

For our purposes it will be convenient to use the origi- 
nal definition, Eq. H3.1|l . The implications of our results 
for technicolor theories would be essentially the same if 
we used the expression H3.2|l . With cither definition, since 
a one-loop heavy fermion correction to the Z propagator 
has a prefactor (g^ -I- (/'^)/(167r^), where g and g' are 
the respective SU(2)i and U(1)y gauge couphngs, and 
since (g^ + g'"^) / (IQir'^) = ae™/(47rs^c^), the prefactor 
is'^c^ /aemi^z) in the definition of S cancels out the 
leading dependence on the electroweak gauge couplings 
(evaluated at the scale mz), yielding a quantity that de- 
pends on the intrinsic properties of the strongly coupled 
SU(iV) gauge theory. 

Now, suppressing the SU{Ntc) gauge index, we write 
the fermions as a vector, ip = {ij^i, ...^ipjsij). We then 
define vector and the axial-vector currents as 



F;(a;) = ^(x)T'^7,,V(x) 
Alix) = ^(x)T^^75V(x) , 



(3.3) 



where the Nf x Nf matrices (a = 1, 7V|- 1) are the 
generators of SU{Nf) with the standard normalization 
Tr(r'^r'') = In terms of these currents, the two- 

point current-current correlation functions flyy and Haa 
are defined via the equations 



8'^' (^-5..)n;.(g^) 



I I d^x e'^- (0|T(J^"(a;)J,^(0))|0) 

(3.4) 

where {x) — {x) , A" (x) . With the above normal- 
ization of T", n(q^) measures the contributions to the 
time-ordered product in Eq. H3.4|l per fermion. Given 
that, to a good approximation, different technifermion 
doublets contribute equally to 5, it is natural to define a 
reduced quantity, S, that represents the contribution to 
S from each such pair, viz.. 



(3.5) 



Then, in terms of the current-current correlation func- 
tions defined above, S, as defined in Eq. (|3.1() . is given 

by 



S = 4tt — [nvv{q^)-nAA{q^)] 



It is convenient to define the compact notation 



n 



V-A 



{q^) = nvv{q^)~TlAA{q^ 



(3.6) 



(3.7) 



As is evident from Eq. (|3.6|) , one may also consider S in 
a different context, namely that of an abstract vectorial 



SU(Af) gauge theory with Nf massless fermions trans- 
forming according to the fundamental representation of 
this group, and with all other interactions much weaker 
in strength than the SU(A^) gauge interaction. In this 
case, in contrast to technicolor models, where Nf is even 
(since the technifermions have left-handed components 
forming SU(2)i doublets), Nf can be even or odd (being 
restricted to be less than Nf^cr so that the theory is in 
the confinement phase with spontaneous chiral symmetry 
breaking). Here, one could naturally define the contribu- 
tion to 5* from each fermion individually, namely, S/Nf. 
However, since our main application will be to techni- 
color, we shall continue, as in Ref. to present our 
results in terms of S. 

Because of the asymptotic freedom of the SU(Af) the- 
ory, for Euclidean much larger than A^, dimensional 
considerations imply that, asymptotically, Hy-Aiq'^) — 
{iptp)'^ / q^ up to logs arising from anomalous dimensions. 
Combining this property with the analytic properties of 
Iiv-A{q^), one can write the following dispersion rela- 
tion, with t = q^: 

t: Jo s-t-ie 
The dimensionless spectral function pj{s) is defined as 



Pjis) 



/m(Hjj(s)) 



ITS 



(3.9) 



This spectral function encodes information about the 
hadronic states produced by the current J. In terms 
of these spectral functions, one then has 



ds 
s 



Pv{s) - Pa{s) 



(3.10) 



An early discussion of the integral on the right-hand side 
of Eq. (|3.10|l and its connection to 7r+ i^i'f j decay 
appeared in |46j| . It is of interest to comment on the 
two equivalent expressions, Eqs. (|3.fi|l and H3.1()|l for S. 
The first of these, Eq. (|3.6|) . expresses S in terms of 
the slope of Ilv-A{q^) at q^ — 0. In contrast, the sec- 
ond, Eq. H3.10|l expresses it as an integral of 1/s times 
the difference of the physical spectral functions for the 
vector and axial vector currents in the timelike region, 
which depends on the spectrum of hadronic states and 
the strengths of their couplings to these currents. Thus, 
naively, one might think that these two expressions de- 
pend on rather different properties of the SU(Af) theory. 
The fact that they are actually equivalent is a conse- 
quence of the analytic properties of Iljj{q^) which are 
used, via the Cauchy theorem, to derive the resultant 
dispersion relation. Thus, although it is evaluated at a 
single point, the derivative ny_^(0) "knows" about the 
full hadronic spectrum, including the presence or absence 
of walking. 

In passing, we note that the spectral functions directly 
in terms of the products of currents, as 



(o|T(y;(x)i^,^(o))|o) 



5 



ds 



and 



ah 



ds 



PA(s){sg^,i, - q^qy) q^q^ 



(3.11) 



2 



(3.12) 



of a*, closer to QCD with Nf = 2 or 3 and found a 
dramatic growth in fp/A, approaching values nearer to 
QCD, as expected. 

Because of the asymptotic freedom and consequent 
asymptotic decay oiIlv-A{q^) for large Euclidean q^ and 
the fact that the fermions are massless in the underlying 
SU(iV) theory, the following integral relations (first and 
second Weinberg sum rules) hold 49]- 52J: 



ds 



Pv{s) - pa{s) 



= \Iv-a{^) = .fp (3.16) 



where one explicitly separates out the contribution from 
the (J = part of) Pa{s) at s = due to the Nambu- 
Goldstone bosons (NGB's). In Eq. p. 121) . the quantity 
Jp (where P stands for "pseudoscalar" ) is the TV^-flavor 
generalization of defined by the transition matrix ele- 
ment 



(0|^^(0)|7r^(<z)) =zg,Jp<5' 



ah 



(3.13) 



with a,b — 1,2,..., Nj — 1. In actual QCD, the chiral 
symmetry is explicitly broken by the u and d current- 
quark masses (and also by electroweak interactions), so 
that the pions decay, and, in particular, the dominant 



weak decay of the 7r+, tt"* 



has a rate propor- 



tional to Thus, fp might be called the generalized 
Nambu-Goldstone decay constant, but we will avoid this 
term, since in our basic SU(iV) theory with other inter- 
actions turned off, these Nambu-Goldstone bosons are 
exactly massless and do not decay. In the chiral limit 
of QCD, with uiu = VTLd = 0, it has been estimated that 
{fTT)ch.iim.l Itt — 0.935 0, so that, with the physical value 
= 92.4 ± 0.3, one infers that {f,r)ch.iim. ^ 86 MeV 
(with a theoretical uncertainty of several per cent from 
the chiral extrapolation). This slight decrease will not be 
important for our work. 

The constant fp may be calculated by first solving the 
Schwinger-Dyson equation for the momentum-dependent 
dynamical fermion mass S(p) and then substituting this 
into the Pagels-Stokar relation [4^ . 



f2 _ ^11 



47r2 



ydy 



[y + mv) ? 



(3.14) 



Calculations using this method |2lL l25l l2a | have shown 
(for N = Nc = i) that as Nf increases from the value 
2 toward Nf^r, the generalized quantity fp de- 



creases, as is expected, since fp is an order parameter for 
spontaneous chiral symmetry breaking. Furthermore, in 



the strong walking limit Nf 



N 



f.nr. I.e., g^, 



has been found that fp vanishes like |3, |4a| 



fp = CfA exp 



\ OLrr / 



-1/2' 



where c/ is a constant. In Ref. {23 we calculated fp in 
the crossover region between this extreme walking limit 
and smaller values of N f^ corresponding to larger values 



and 



ds s 



pvis) - pa{s) 







(3.17) 



Clearly, Eqs. IpTTnjl . (pTTHjl . and IpTTTjl are of similar 
form, viz., integrals of pvis) — Pa{s) weighted, respec- 
tively by s~^, 1, and s. Because of these different factors, 
the contributions to the integral (|3.1U|) for S are weighted 
toward smaller s values, while the contributions to the 
integral H3.17|l are weighted toward larger s values, as 
compared with those for the integral in Eq. (|3.16() . 

Finally, in the context of an abstract SU(iV) gauge 
theory (with other interactions turned off), it is worth- 
while to mention the connection between S and the coef- 
ficient denoted Liq. To explain this connection, we recall 
the definition of Liq. Provided that Nf > 2, the the- 
ory contains N'j — 1 massless Nambu-Goldstone bosons. 
(For this discussion, we turn off electroweak interactions 
completely; for our application to technicolor, we do not 
turn them off, and hence three of the would-be Nambu- 
Goldstone bosons become the longitudinal modes of the 
and Z. Furthermore, in the TC/ETC context, other 
would-be NGB's gain masses from color and ETC inter- 
actions that explicitly break the SU(7V/)l x S\]{Nf)p 
global chiral symmetry.) A useful description of the low- 
energy physics for energies E <^ ^''[Ip is then provided 
by a chiral Lagrangian (see. e.g., Q and references to 
earlier work therein). This is a function of the chiral 
fields 



U = exp 



N —^ 



(3.18) 



For example, for QCD with just the two light quarks u 
and d, U would have the form U ~ e^'^-'^/f^. Now let 
us denote the elements of the global symmetry groups 
SU(7V/)l and SV{Nf)p as Ul and Up. Then under a 
transformation by an element of the chiral symmetry 
group SU(iV/)L X SV{Nf)R, U UlUUI- One can also 
formally introduce external gauge fields (contracted with 

(3.15) 

the respective generators of SlJ{Nf)L and SU(7V/)_r), 



and R^. These have the transformation properties Lf^ 



ULL^Ul-i{d^UL)Ul and ^ URR^Uj^-iid^URpl 
In terms of these, one constructs the covariant derivative 



iURu- One also defines external 
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field strengths {Fl)^^, and (F^)^^ in the usual manner. 
The lowest-order term in this effective Lagrangian is then 



Tv[{D^U)Hd^'U)\ 



(3.19) 



In terms of these quantities, L-^q is the coefficient of the 
term 



Oio -Tr(C/t(Fi)^,C/(Fflr), 



(3.20) 



where the superscript (r) refers to the renormalized quan- 
tity. From this, one obtains a quantity which is con- 
structed to be independent of the renormalization scale 



1 



In 



The relation with S is then 



5 



-167rLi 



(3.21) 



(3.22) 



Note that although Lio requires Nf > 2 to be well- 
defined, it does not require Nf to be even. For the specific 
case of QCD, fits to experimental data, including, in par- 
ticular, the radiative decay 7r+ e^u^j, yield the value 



5 = 0.33 ±0.04 



(3.23) 



To the extent that this is dominated by the two light 
quarks u and d, one has Nf = 2, and hence N^ = 1, so 
that the measured value of S for QCD is also equal to 
S. The fact that the light-quark vector and axial-vector 
mesons p and oi largely saturate the expression for S, 
Eq. (|3.10|l . is consistent with this conclusion. An ap- 
proximate calculation of S has been performed using the 
ladder approximation to the Schwinger-Dyson and (inho- 
mogeneous) Bethe-Salpeter equations for QCD {N = 3) 
with Nf — 2 quarks of negligible mass . Studies have 
also been done for the case where one neglects the strange 
quark mass m^, i.e., N = 3, Nf = 3 0,13. Since for 
either of these values of Nf the beta function of the QCD 
theory does not exhibit an infrared fixed point, it is nec- 
essary to cut off the growth of the strong coupling. For 
typical cutoffs, it was found that the SD-IBS calculations 
tended to yield slightly too large a value of S' = 5', namely 
S ~ 0.45 — 0.5 ,21^^], rather than a value in the la ex- 
perimental range 0.29 ^ S ^ 0.37. This suggests that 
in the SD-IBS approach, used for a vectorial confining 
SU(iV) gauge theory with small values of Nf such that 
the theory has no perturbative IR fixed point, with a 
typical IR cutoff on the coupling, tends to overestimate 
S by about a factor of 1.4. Since the calculation of S 
in QCD is a problem in strongly coupled, nonperturba- 
tive physics, and the calculational method that was used 
is only approximate (neglecting, for example, instanton 
contributions), one should probably not be surprised that 
it does not precisely reproduce the measured value of S. 



IV. CURRENT-CURRENT CORRELATION 
FUNCTIONS IN TERMS OF BETHE-SALPETER 
AMPLITUDES 

In this section, we explain how the current-current cor- 
relation functions are obtained from the Bethe-Salpeter 
amplitudes, which will be calculated via the inhomo- 
geneous Bethe-Salpeter (IBS) equation 53J- 63], [48j| . 
These Bethe-Salpeter amplitudes x^'^\ where J = V or 
A, are essentially form factors, whose behavior in the 
timelike region describes the coupling of the given cur- 
rent to physical hadronic bound states that can be pro- 
duced by this current, with an analytic continuation into 
the spacelike region. Here we will only need these am- 
plitudes for the spacelike region < and at the point 
= 0. The amplitudes may be defined in terms of the 
three-point vertex function as 



(27r) 

= 6^1 d^xe^'^-^0\T{i^'^f{r/2) i^,f>p{~r/2) J«(a:))|0), 

(4.1) 

where {f,f'), {j,k) and {a, (3) are, respectively, the fla- 
vor, gauge, and spinor indices. Closing the fermion legs 
of the above three-point vertex function and taking the 
limit r 0, one can express the current correlator in 
terms of the Bethe-Salpeter amplitude as 




(e.G('^))x(-^)(p;<Z,e) 



where \64 



7m > 



7m 75, 



(4.2) 



(4.3) 



and an average has been taken over the polarizations, so 
that Iljj{q'^) does not depend on the polarization e. 

We expand the Bethe-Salpeter amplitude xi^''(p;9j£) 
in terms of a complete bispinor basis with basis elements 
r^'^'' and the invariant amplitudes xi'^^ 



where we define 



Q/3 



— L J ap 



i=l 



so that > in the spacelike region, and set 



9m 



_ 9m 



(4.4) 
(4.5) 

(4.6) 
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The bispinor basis elements can be chosen in such a man- 
ner that they have the same spin, parity and charge con- 
jugation as the corresponding current J^{x). For the 
vector vertex we adopt the foUowing bispinor basis ele- 
ments: 



1 



(V) _ 



1, 



^ 3! 



TP^i{p-q){e-ph TP = ^[^,i]{e-p), (4.7) 

where [a,6, c] = a[fe, c] + 6[c, a] + c[a, 6]. For the axial- 
vector vertex we use the bispinor basis elements 



ri^'^^-i^JKe-p) ip-q) 75. 



5 



i){e ■ p) 75, T\ ' ^ ^ {e-p) [p- q) 75, 



(4.8) 



Given the charge conjugation properties of the vector and 
axial- vector currents and the above choice of the bispinor 
basis elements, it follows that invariant amplitudes xP 
are even functions of p • q. 

In the present analysis it is convenient to choose the 
Lorentz reference frame so that only the timelike compo- 
nent of is nonzero. Since we are working in the space- 
like region (which avoids physical mass singularities), we 
thus use a Wick rotation with 



g^ = (*Q, 0,0,0). 



(4.9) 



Similarly, for the relative momentum p^ we perform a 
Wick rotation and parametrize it in terms of the real 
variables u and w (with dimensions of mass) as 



m 2 2 2 

p ■ q ^ —Q u , p ^ —u — w . 



(4.10) 



Hence, the invariant amplitudes x- are functions of u 
and w: 



Xi = Xi [u,w;Q). 



(4.11) 



Owing to the charge-conjugation properties of the Bethe- 
Salpeter amplitude x^"^-* and the bispinor basis elements 
defined above, the invariant amplitudes x^Pi''^^^) satis- 
fies the relation 



Xy\u,w;Q) = xf'i-'^^'^'-'Q) 



(J) I 



(4.12) 



Using this property of the invariant amplitudes, we 
rewrite Eq. H4.2|l as 



N 1'°° r 

nvy(Q') = ^ / du 

Jo Jo 



dw 



^xP{u,w;Q) + ^xP{u,w;Q) 



(4.13) 



nAA(Q') = ^ du 

^ Ja Jo 



dw 



xi {u,w\Q)~- —Xq '{u,w;Q) 



(4.14) 



Here, we have used the expanded form of the Bethe- 
Salpeter amplitude given in Eq. (|4.4|l and carried out 
the three-dimensional angular integration. 
From Eqs. ifirr^ and (jlT^ . it follows that 



n 



V-A 



1 /N 
3 V'2" 



E 



d'^p 



Tr 



i(27r)4 

^x''''\p;q,e)- ^j5X^^^{p;q,e) 



— i du 



dw w'^ 



xP {u, w; Q) + xP {u, w; Q) 



w 



+^\xP\u,w-Q)+xnu^w-Q) 



(4.15) 

Ahhough both Iivv{Q'^) and Haa(Q^) are individu- 
ally logarithmically divergent, the underlying 'S\]{Nf)L x 
S\]{Nf)ii chiral symmetry guarantees that these diver- 
gences cancel in the difference Hyy — n^i^, which is there- 
fore finite. 



V. INHOMOGENEOUS BETHE-SALPETER 
EQUATION 

In this section we discuss the full (inhomogeneous) 
Bethe-Salpeter equation, which we will use to calculate 
current-current correlation functions and, from these, S 
(and, as a check, also fp). The IBS equation is a self- 
consistent description of the coupling of a current J° 
to fermion-antifermion bound states. This coupling is 
represented by the Bethe-Salpeter amplitude x^''^- The 
four-momenta assigned to the fermion and antifermion 
are p^ = p + {q/2) and p^ = —p + {q/2) so that the to- 
tal momentum of the bound state is q, and the relative 
momentum of the ip and ijj is 2p. Since we are deal- 
ing with J = 1 bound states, the bound-state amplitude 
also depends on the polarization vector e, which satisfies 
e ■ q = and e • e = — 1. A grapical indication of the in- 
homogeneous Bethe-Salpeter equation structure is given 
in Fig. n The IBS equation is 



T{p;q) x'^''Hp;q,£) 



■G(^) + K[p-k)*x''^\k-q,€). 

(5.1) 
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VliP-.'Dt 




P+I/2 



+ 




p-y2 



P-% k-1/2 



FIG. 1: A graphical expression of the IBS equation in the 
(improved) ladder approximation. 



Here the kinetic part T is 

T{p- q) = -Sp\p + q/2) ® Sp\p - q/2) , (5.2) 

where Spip) = 1/ {A{p)^—T,{p)) . We follow the standard 
procedure of using the Landau gauge in calculations with 
the Schwinger-Dyson and Bethe-Salpeter equations since 
in this gauge the fcrmion wave function renormalization 
factor A{p) = 1. The physical results are, of course, 
gauge- invariant. The Bethc-Salpeter kernel K in the im- 
proved ladder approximation is expressed as 



Kip; k) 



2/ 



{p~kf 

{p-k)^{p~k)^ 



[p - kf 



•7''®7^ (5.3) 



where C2/ was given in Eq. H2.5|) . In the above expres- 
sions we use the tensor product notation 



and the inner product notation 

cf^k 



(5.4) 



K{p;k)*x^'\k;q,e)^ 



i(27r) 



K{p,k)x{k;q) ■ (5.5) 



where summations over Dirac indices are understood. (In 
contrast, for our previous calculations 25, 26J of meson 
masses, we only needed to use the homogeneous Bethe- 
Salpeter equation.) 

As in Refs. 24, EllH IHiOl , we make the ansatz for 
the running coupling, after Euclidean rotation. 



aiPE,qE) 



(5.6) 



where the subscript denotes Euclidean. Since a would 
naturally depend on the giuon momentum squared, (p — 
q)^ = p^ +q^ — 2p-q, the functional form H5.6() amounts to 
dropping the scalar product term, —2p ■ q. As discussed 
in Ref . |25j | , this is a particularly reasonable approxima- 
tion in the case of a walking gauge theory because most 
of the contribution to the integral comes from a region of 
Euclidean momenta where a is nearly constant. Hence, 
the shift upward or downward due to the —2p ■ q term 
in the argument of a has very little effect on the value 
of this coupling for the range of momenta that make the 
most important contribution to the integral. The ap- 
proximation H5.6|l enables one to carry out the angular 
integration analytically. 



nv_A(Q')//p 



1 1 


1 1 1 

Otil! = 1 .0 




a* = 1.2 




a* = 1.4 




a* = 1.6 




-.^^^ a* = 1.8 




. actual QCD 


1 1 





0.2 0.4 



qV/p^ 



FIG. 2: Plotofnv-A(0^)//p as a function of in the space- 
like interval < QV/p < 1-2, for a* = 1.0, 1.2, 1.4, 1.6, 1.8. 
As indicated, the horizontal axis refers to the quantity / fp. 
For comparison, Hy-AiQ^)/ fp for QCD with Nf = 3 mass- 
less quarks is also plotted. See text for further details. 



The momentum-dependent dynamical mass S(p) for 
the fermion is obtained from the Schwinger-Dyson equa- 
tion. 



i:{p) = -K{p,k)*SF{p). 



(5.7) 



We use the same kernel K{p, k) here as in the IBS equa- 
tion in order to respect the SU(7V/)l x SU(iV/)i?^ chiral 
symmetry [g^, |^ • The numerical method that is used 
for solving the SD equation and the IBS equation involves 
approximating the respective integrals by discrete sums 
and is the same as in Ref. HH ■ The reader is referred 
to these papers for more details on this method. 



VI. RESULTS AND DISCUSSION 

In this section we present the results of our calculations 
of I{v-a{Q^) and S using the Schwinger-Dyson and (in- 
homogeneous) Bethe-Salpeter equations. It is appropri- 
ate to include an obvious cautionary remark that these 
calculations involve strong couplings a of order unity, and 
therefore there could be significant corrections to the (im- 
proved) ladder approximation used in our solutions of the 
Schwinger-Dyson and Bethe-Salpeter equations. Accord- 
ingly, as one check on the reliability of our methods, we 
have also carried out a comparison of fp obtained from 
the SD and IBS equations via eq. (|3.16l) with fp obtained 
from the SD equation via the Pagels-Stokar relation. 

In Fig[51we plot our calculated values oiWy - a{Q^) I fp 



for 



1.8, 1.6, 1.4, 1.2 and 1.0, as functions of the di- 



mensionless quantity Q"^ / f%- The slope of each curve at 
= is equal to ~S/ {Att) for the given value of a* . We 
first observe that I{v-a{Q'^) is almost linear at = 0, 
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A 

Sn 



ttf 



» o o o 



. o <> . . . . o o 



1.6 1.8 a * 



a * 



FIG. 3: Plot of Sn for several values of a, in the range 0.9 < 
a* < 1.8. As indicated by the subscript n, the values are 
normalized by the value of 5 at a* = 1.8. 



FIG. 4: Plot of the ratio {fp)ps/(fp)wi, where [fp)ps and 
{fp)wi denote fp calculated via Eq. !);-{. 14|l and via Eq. 
l|3.16|l . respectively, as a function of a,. 



with a small positive second derivative. This is the justi- 
fication for our earlier statement that the implications of 
our findings for technicolor would be essentially the same 
whether we used the definition of S in terms of a deriva- 
tive at = or a finite difference, analogously to Eq. 
(|3.2(l . Secondly, from Fig. |21it is clear that the magni- 
tude of the slope, and hence S, decreases as a* decreases 
toward the chiral phase transition point acr- 

In Fig. 131 we plot the values of S, derived from the 
slope of IIv-a{Q^) at — 0, for several values of a*. 
As indicated by the subscript n, the values are normal- 
ized by the value of 5 at a* = 1.8. This figure shows that 
S decreases by about 40 % as a* is reduced from 1.8 to 
0.9, or equivalently (c.f. Eq. (|2.3|l 'l as Nf is increased 
from 10.3 to 11.6. As our calculation of meson masses 
in Ref. showed, this is a crossover region, in which 
the theory is changing from QCD-like, non-walking be- 
havior at smaller Nf to the walking regime at larger Nf 
approaching Nf ^r- Reinserting the factor of N]j = Nf/2 
to get S itself, we obtain a decrease by about 30 % in S, 
since No only increases by about 10 % over this range. 
Thus, our calculation shows that for this range of values, 
S decreases significantly as one moves from the QCD- 
like to the walking regimes. This finding is an important 
result of our present study. 

As a check on our calculational methods, we com- 
pare fp calculated in two different ways: via the Pagels- 
Stokar relation, Eq. (|3.14|l . denoted {fp)ps, and by 
the first Weinberg sum rule (Wl) or equivalent relation 
Uv-AiO) = fp in Eq. denoted {fp)wi- In view 

of the fact that the Pagels-Stokar itself is approximate 
and that our solution of the SD and IBS equations in- 
volves the ladder approximation and the neglect of com- 
pletely nonperturbative contributions such as those due 
to instantons, we do not expect exact agreement between 



these two different methods of calculation. In Fig. 0]we 
present a plot of the ratio {fp)ps I {fp)wi as a function 
of a* . The closeness of this ratio to unity gives one mea- 
sure of the accuracy and reliability of our calculations. 
We see that the ratio is within about 20 % of unity and 
is essentially independent of a* in the range considered, 
with the Pagels-Stokar method yielding a slightly smaller 
value than the expression in terms of the current-current 
correlation functions. This gives us further confidence in 
the results of our SD-IBS calculation of S. 

Our finding that S and S are reduced as one moves 
from the QCD-like regime toward the walking regime 
of an SU(A^) gauge theory is in agreement with the ap- 
proximate analytic results of Refs. I20L |22L I23il |. and 
it complements those works, being based on a numeri- 
cal solution of the Schwinger-Dyson and Bethe-Salpeter 
equations. In the sub-interval 0.9 < a* < 1.0 closer 
to the walking limit our results coincide with those in 
Ref. [2^. Our use of a larger interval has the advan- 
tage that we are able to observe a larger reduction in 
S as a* decreases than was done in Ref. 2M- We have 
not attempted here to examine the extreme walking limit 
(a* — acr)lcy.cr ^0"*". As discussed in Ref. 0, it be- 
comes increasingly difficult to obtain an accurate numer- 
ical solution for S in this limit because of very strong 
cancellations. 

In addition to this decreasing trend of S, one may 
also discuss the absolute magnitude of S. Our calcu- 
lation yields S = 0.47 at a* = 1.8. If S continues to 
be a monotonically function of Nf (and hence, in the 
range of interest here, also a monotonic function of a*), 
then an extrapolation of our calculation to the QCD case 
with N f = 2 or, with the strange quark mass neglected, 
Nf = 3, would predict a value of 5' ^ 0.5. This is similar 
to the value that was obtained in earlier studies using 
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the SD and IBS equations in a different manner than 
here, where it was necessary to introduce an cutoff on 
the growth of the strong couphng in the infrared [2lll2^ . 
In Ref. ^2§\ it was shown that if one used a cutoff that 
led to a very large value of the coupling, one could get a 
result for S in agreement with the experimental value 
(|3.23(l . but the reliability of the calculational method 
in the presence of such a large coupling was not clear. 
We shall adopt the optimisitc viewpoint here of giving 
greater weight to the change in 5 as a function of Nf than 
to the absolute value of S itself. Equivalcntly, one could 
envision applying an overall correction factor of about 
2/3 to the absolute magnitude of S so that the value for 
small Nf matches that in QCD. Physically, this factor 
would be regarded as correcting for the strong-coupling 
effects not included in the SD-IBS analysis. 

Although one cannot use perturbation theory reli- 
ably to calculate S* in a strongly coupled gauge the- 
ory, the perturbative result is often used in discussions 
of constraints on new physics, and hence it is worth- 
while to see how our results compare with the pertur- 
bative computation. A one-loop perturbative calcula- 
tion with degenerate fermions having effective masses 
satisfying {2'E /mz)'^ ^ 1 yields the well-known result 
Spert. = ND,tot./{QT') where here ND,tot. = NdN, i.e, 



5. 



pert. 



N 



(6.1) 



In QCD with Nf ^2 and N = Nc^i, this perturbative 
calculation would predict SqcD.pert. — l/(27r) ~ 0.16. 
The experimental value in Eq. (|3.23|) is approximately 
twice this; Sqcd — '^SqcD.pert.- The reductions that we 
have found in S and S for the range of a* investigated 
suggest that in a walking theory, much or all of the above 
factor of 2 might be removed, and the true value of S 
might well be comparable to, or, indeed, perhaps less 
than, the perturbative estimate. 

For reference, in the one-doublet and one-family tech- 
nicolor models, the perturbative expressions for the tech- 
niparticle contributions to S are Srcpert. = A^Tc/(67r) 
and Srcpert. = 2iVrc'/(37r), respectively. With Ntc = 
2, these take the values SQCD,pert. — 0.1 and 0.4, respec- 
tively. Fits to precision electroweak data yield allowed 
regions in S and two other parameters describing mod- 
ifications of the Z and W propagators by new physics 
beyond the Standard Model, namely the parameter T 
measuring violations of custodial SU(2) from this new 
physics and a third parameter, U, of somewhat less im- 
portance here. Since the Standard Model expression for 
S includes a term {l/{6'K))lii(mH /rriH^ref.), the resul- 
tant allowed regions depend on the choice of the reference 
value of the SM Higgs mass, mn^ref.- The comparison of 
these with a technicolor theory is complicated by the fact 
that technicolor has no fundamental Higgs field; some- 
times one formally uses rnH,ref. ~ 1 TeV for a rough 
estimate, since the SM with tuh ^ 1 TeV has strong 
longitudinal vector boson scattering, as does technicolor. 
However, this may involve some double-counting when 



one also includes contributions to S from technifermions, 
whose interactions and bound states (e.g., techni- vector 
mesons) are responsible for the strong W^W£ and Z^Zl 
scattering in a technicolor framework. The current fit 
[TgI It^ disfavors values of 5 ^ 0.2. Our findings in 
this paper suggest, in agreement with the previous works 
noted above using different methods, that the constraint 
on walking technicolor models could be less severe than 
would be inferred from the perturbative formula for the 
technifcrmion contribution to S. 



VII. SUMMARY 

In summary, using numerical solutions of the 
Schwinger-Dyson and (inhomogeneous) Bethe-Salpeter 
equation, we have calculated the S parameter as a func- 
tion of the approximate infrared fixed point, a*, or equiv- 
alently, the number of massless fermions, Nf, in a vec- 
torial, confining SU(Af) gauge theory. We have focused 
on the crossover region between the walking and QCD- 
like (non-walking) regimes. Our results show that S and 
also S decrease significantly as a* decreases in this range. 
This trend agrees with earlier indications of a decrease 
in S in walking gauge theories. We have discussed the 
implications for technicolor models. 
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thanks T. Appelquist for useful comments. This re- 
search was partially supported by the grant NSF-PHY- 
03-54776. 



VIII. APPENDIX 

In this appendix we present some results on vector 
meson dominance (VMD) and their relevance to the 
crossover region that we are studying numerically. In 
QCD itself, VMD has served as a useful approximate 
model. The matrix elements for the production of a vec- 
tor or axial-vector meson from the vacuum are 



{p^\V'\0)^^6^'f, 



{a1\Al\0)=tS^''fa,e, 



(8.1) 



■2) 



where are the respective polarization four- vectors. In 
QCD, fp ~ 154 MeV, so that the first Weinberg sum rule 



yields = - /2 = 123 MeV. A pole-dominated 

form for the current-current correlation functions in the 
(simplistic) narrow-resonance approximation is 



my — — it m\ — — it 



(8.3) 



where, for QCD, the ground-state vector and axial- vector 
mesons are the p(776) and ai(1230), i.e., V = p and 
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A = ai. Using the principal value relation 

= P (^—] ± iTTS{s - So) , (8.4) 

the pole-dominated form (|8.3|) gives, for the spectral 
functions, 

Pvis) = f^S{s - ir?y) (8.5) 

and 

p^(s) = /l<5(s-mi) (8.6) 

Although these are useful simplifications, more realistic 
analyses of the spectral function sum rules take into ac- 
count the finite widths of the p and a\ resonances, and 
also nonresonant contributions 52]. 

Expanding Eq. (|8.3() for large Euclidean g^, one ob- 
tains 

n._.(,-)^^ ^^^<;/>^) +o(^l). (8.7) 

The asymptotic freedom of the SU(iV) gauge theory im- 
plies that for \q^\ > A^, nv_^(q^) cx l/g'*, so this yields 
the relation 

fWv - !Wa = , (8.8) 

which is the second Weinberg sum rule, evaluated in 
the VMD approximation. Substituting the forms for the 
spectral functions into Eq. H3.10|) yields 

SvMD = -167rLio = 47r ( ^ - i|- I (8.9) 

Similar substitutions into the first and second Weinberg 
sum rules yield the relations 

fv-fA = fp (8.10) 

and Eq. H8.8|l . The narrow- width VMD relations are 
approximate, since they neglect the sizeable widths of 
the p and ai , but the Weinberg relations follow from the 
asymptotic freedom of QCD when one sets vriu = md = 0. 

As Nf increases and there is a crossover in the SU(A^) 
theory from QCD-like behavior to walking behavior, 
there are resultant changes in how these integral rela- 
tions are satisfied. For both the QCD-like and walking 
regimes, the underlying asymptotic freedom of the the- 
ory implies that the current-current correlation functions 



Iiv-A{q^) have a asymptotic falloff (up to logs) for 
Euclidean 3> A^. However, in a walking regime, the 
scale of chiral symmetry breaking, given by /p ^ S, is 
much less than A, and in the interval <^ <^ h? , 
liv-Aiff') have a less rapid falloff, ~ l/g^^^'*'), where 
7 is the anomalous dimension of the bilinear fermion 
operator and is 7 ~ 1 for strong walking, whence a 
njj(9^) 1/9^ falloff in this interval. Physically, the 
coupling a is strong, of order 0(1) but slowly varying 
for an extended interval of energies between E and A. 
Thus it is plausible that in a walking theory, rather than 
just the lowest-lying vector and axial-vector resonances 
making important contributions to the various spectral 
function integrals, there could be significant contribu- 
tions from a number of higher-lying states with the same 
quantum numbers [SOj . In this case, still retaining the 
VDM approximation, one would generalize the expres- 
sions above to be sums over these higher-lying states. 
For example, the evaluation of Eqs. (|3.10|l . H3.16|l . and 
(|3.17|) in the VMD approximation would read 

5,,„.4,£(|--i|-) (8,11) 
= (8-12) 

i 

and 

E(/v.'™y.-/l'^A.) = 0. (8.13) 

i 

More generally, however, it is not clear that the VMD 
model would apply in this walking regime. One analytic 
approximation is to use VMD for a partial evaluation 
representing contributions from ^/s ~ S and supplement 
this with a term due to a fermion loop j22l | . In the region 
that we have concentrated on in this paper, in which the 
theory is crossing over between strong walking behavior 
and QCD-like behavior, one expects that the saturation 
of the spectral function integrals has a form that is in- 
termediate between walking and QCD-like. In our actual 
calculation of S via Eq. H3.6|l . we only need to calculate 
ny_^(0), which we do by means of numerical solutions 
of the relevant SD and IBS equations, so we do not have 
to deal with questions of how the spectral functions be- 
have in a walking gauge theory. In future work it would 
be worthwhile to use the connection between these two 
different ways of calculating S to undersand the hadronic 
spectrum better in such a walking gauge theory. 
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